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We present a novel mechanism for thermalizing a system of particles in equilibrium and nonequi- 
librium situations, based on specifically modeling energy transfer at the boundaries via a microscopic 
collision process. We apply our method to the periodic Lorentz gas, where a point particle moves 
diffusively through an ensemble of hard disks arranged on a triangular lattice. First, collision rules 
are defined for this system in thermal equilibrium. They determine the velocity of the moving 
particle such that the system is deterministic, time reversible, and microcanonical. These collision 
rules can systematically be adapted to the case where one associates arbitrarily many degrees of 
freedom to the disk, which here acts as a boundary. Subsequently, the system is investigated in 
nonequilibrium situations by applying an external field. We show that in the limit where the disk is 
&X), endowed by infinitely many degrees of freedom it acts as a thermal reservoir yielding a well-defined 

^ ■ nonequilibrium steady state. The characteristic properties of this state, as obtained from computer 

' simulations, are finally compared to the ones of the so-called Gaussian thermostated driven Lorentz 

^ : gas. 



> 



On 
0\ 



I 



X 



PACS numbers; Q5.20.-y, 05.45.+b, 05.6Q.+W, 05.7Q.Ln, 44.90.+C 



I. INTRODUCTION 
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, The investigation of transport properties of many-particle systems in nonequilibrium situations generally requires 
thermostats which remove excess energy to ensure the existence of nonequilibrium steady states with constant, or on 
average constant, energy Hoover, Evans, Nose and others developed methods of thermostating by introducing a 

Ch ■ momentum-dependent friction coefRcient into the microscopic equations of motion, modeling the interaction of parti- 
cles with a thermal reservoir P,p|J^~[l^]. These methods are deterministic and time reversible, in contrast to stochastic 
thermostats [p|,p^ [T^] . In this paper we propose and analyze in detail an alternative deterministic thermostat, based 
Q on including energy transfer for a microscopic collision process between the moving particles and the boundaries 
instead of using a momentum-dependent friction coefficient. A short account of the main idea was reported in [ p7| . 

The two basic versions of a conventional deterministic thermostat are the Gaussian thermostat and the Nose- 
Hoover thermostat. The Gaussian thermostat |^|-|lo| creates a microcanonical ensemble in equilibrium and keeps the 
total energy (isoenergetic), or the kinetic energy (isokinetic), constant in nonequilibrium. The equations of motion of 
this thermostat can be derived from Gauss' principle of least constraint [tIJioI]. The Nose-Hoover thermostat 
creates a canonical ensemble in thermal equilibrium and keeps the energy on average constant in nonequilibrium. 

Though the microscopic equations of deterministic thermostated systems are time reversible the macroscopic dynam- 
ics is irreversible in nonequilibrium leading to momentum and energy fluxes with well-defined transport coefficients 
|p|,P|JT^|20|] . Macroscopic irreversibility is based on the fact that only one direction in time is dynamically stable in 
these systems, whereas the time reversed direction is dynamically unstable [Eo|,pl|. This implies a contraction of 
the phase space onto a fractal attractor during the forward evolution | |l^ , p2| -p7|P In contrast, purely stochastic ther- 
mostats are expected to typically lead to a smooth phase space density in nonequilibrium ||l6|,^. In agreement with 
the phase space contraction, it has been found that the sum of the Lyapunov exponents is negative in deterministic 
thermostated systems. It has been argued by many authors that the rate of this phase space contraction is related 
to the thermodynamic entropy production p0| , p9| -^ . Furthermore, relations between the sum of the Lyapunov ex- 
ponents and the corresponding transport coefficient have been derived |l9|j30| , |3^ , ^ , ^ -|4^ . In short, deterministic 
thermostats provide an important approach to modeling nonequilibrium steady states and establishing interesting 
links between dynamical system theory and statistical mechanics PJs],^^^ . 

On the other hand, conventional thermostats are based on a drastic modification of the microscopic equations 
of motion by including momentum-dependent friction coefficients, which implies that the microscopic equations 
cannot be Hamiltonian anymore in their usual physical coordinates. Although there exist methods to relate them 
to generalized Hamiltonian systems by noncanonical transformations [|7|,^-pl| the question still remains whether the 
results obtained from these deterministic thermostats provide general characteristics of nonequilibrium steady states, 
or whether they depend on this particular way of thermostating pq-M ■ 
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As an alternative, a specific mechanism to simulate a steady shear flow without using a thermostat of the above 
kind has been studied in Refs. [|34y37|j5^ . Here, the collision of a particle with a wall is described by rules which change 
the scattering angle but not the absolute value of the velocity of the particle. Open systems with fixed concentration 
gradients at the boundaries have been the subject of another approach to create a nonequilibrium steady state [^3|-^6|. 
A possible link of this approach to thermostated systems is discussed in Refs. ]35|j3^ . 

A simple deterministic one-particle system in which there is evident need of thermostating is the field driven periodic 
Lorentz gas. We recall that the original Lorentz gas model consists of a system of randomly distributed hard disks and 
a particle that moves freely between successive elastic collisions with the disks [ST) . Later, a periodic configuration 
of disks onto a triangular lattice known as the periodic Lorentz gas has been introduced, and still serves as 
a standard model in the field of chaos and transport [^ , ^6|j5^ -|65|| . In case of the driven Lorentz gas, an external 
electric field drives the system into nonequlibrium by accelerating the moving particle while pumping at the same 
time energy into the system through Joule heating. A number of authors developed mechanisms for removing energy 
from this system through a Gaussian isokinetic thermostat, which creates a nonequilibrium steady state with constant 



nonequilibrium as other many-particle systems: the phase space density contracts onto a fractal attractor ||19|,|66|-|68| , 
and the sum of Lyapunov exponents is negative. Relations were derived between this sum of Lyapunov exponents, 
the conductivity and the irreversible entropy production of this system. ||3^ , |33| , ^ - ^ , ^7||6S[| . We note that a model 
almost identical to the driven periodic Lorentz gas, except for some geometric restrictions, is the Galton board, which 
has been invented in 1873 to study probability distributions [ p6| . 

In the present paper we introduce an alternative method of deterministic thermostating which is free of addition 
of new terms in the equations of motion, and illustrate it on the periodic Lorentz gas. The paper is organized as 
follows. In Section II we introduce our model in equilibrium. We define collision rules for the particle which change 
its velocity at a collision with a scatterer such that the dynamics of the system is deterministic, time reversible, and 
yields the microcanonical density in equilibrium. First the disk is equipped with one degree of freedom, and energy 
conservation between particle and disk determines the amount of energy on the disk after a collision. Then the model 
is modified by pretending that the disk has arbitrarily many degrees of freedom. This involves a redefinition of the 
microscopic scattering rules to keep the dynamics microcanonical. The reduced densities of an arbitrarily dimensional 
microcanonical system are then calculated. In Section III we numerically investigate the system in nonequilibrium 
by switching on the external field. We show that in the limit of associating infinitely many degrees of freedom 
to the disk our mechanism keeps the energy of the moving particle constant on average, thus leading to a well- 
defined nonequilibrium steady state. This confirms that our mechanism yields indeed a proper thermostating. The 
characteristic features of the resulting nonequilibrium steady state for our model are discussed explicitly, especially in 
comparison to the Gaussian thermostated periodic Lorentz gas. A summary with main conclusions is given in Section 



Owing to the periodicity of the lattice, it will be sufficient to study the dynamics in one Lorentz gas cell with 
periodic boundary conditions, see Fig. |^(a). As the radius of the disk we take r = 1. For the spacing between two 
neighboring disks we choose w ~ 0.2361, as is standard in the literature to ensure that a diffusion coefficient exists 



The variables intervening in the dynamics are defined in Fig. |l](b): (3 is the angular coordinate of the point 

at which the particle collides with the disk, 7 is the angle of incidence at this point, and a is the angle of flight of the 
particle. The particle has two velocity components, v = (vx,Vy). To establish the energy transfer between particle 
and disk we also assign to the disk a velocity k. Thus, in total our model has three degrees of freedom. To proceed 
further we now need to introduce specific collision rules for the moving particle, which map v and 7 onto v' and 7'. 
Energy conservation between particle and disk then yields the velocity k' onto the disk after a collision. 



We first present a model based on a very simple physical mechanism for a possible energy transfer between particle 
and disk at a collision. This model has the required properties only in a limited parameter range, but it contains some 
basic ideas for the formulation of our thermostating mechanism, which we will introduce in full detail afterwards. 

The position of the disk is held fixed. The velocity of the moving particle at a collision can be split into a normal 
component w„ and into a tangent component Vt- We now interpret the velocity k of the disk as a rotational degree 




IV. 



II. THE MODEL AND ITS EQUILIBRIUM PROPERTIES 




A. Rotating disk model 
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of freedom. This allows us to define a transfer of kinetic energy between vt and k based straightforwardly on energy 
and momentum conservation, while u„ is elastically reflected. The collision rules thus read 

, {m - mk)vt + 2mkk 

m + TUk 

, (irik — m)k + 2mvt 

k ~ 

m + m^; 

v'n = -Vn , 

where m is the mass of the particle and the mass of the disk. Computer simulations show that this model yields a 
microcanonical probability density for a total kinetic energy oi E — 0.5, m = ruk — 1 and large spacings between two 
neighboring disks w > 2.0. However, by decreasing w the probability density starts to deviate from the microcanonical 
density because the particle gets more and more trapped in parts of the Lorentz gas cell. By increasing the mass of 
the disk ruk deviations from the microcanonical density appear already for larger w. 

Thus, only for a certain choice of parameters is this simple approach leading to the desired result, which is that 
the system is microcanonical and shows equipartitioning of energy in all degrees of freedom, while, in general, the 
dynamics is apparently more complicated. In the following we want to investigate whether by introducing more 
generic collision rules we can achieve that the dynamics is microcanonical for any choice of respective parameters. 
Specifically the fact that the energy and momentum conservation law is a linear two-dimensional map of the form 
(uj, k') = f{vt, k) motivates us to define the collision process by a simple, chaotic two-dimensional map as discussed 
amply in the next section. More details of the rotating disk model will be reported elsewhere Jzof . 

B. Modeling the collision process by a two dimensional map 

1. The baker map 

We choose the well-known baker map ||4^ , which we apply to the variables {xb, Vb) — (sin I7I, v). Here Xf, = sin(|7|) 
is the Birkhoff coordinate]^ of 7 leading to £i(sin(|7)|) = 1 |^^. The change of these variables at a collision to 
(xb, Vb') = (sin 17'!,^') (see Fig. |) is thus given by 

{Xb, Vb) = M{xb, Vb) = I (2a;, - l,(y, + l)/2) , Xb > 0.5 ' 

where < < 1. k' is then obtained from energy conservation. Since k is not explicitly contained anymore in 
the collision rules given by Eq. (|l|), one can argue that the detailed dynamics of k is no longer relevant for the 
moving particle. In particular, k need no more be associated to a rotational degree of freedom. Based on our general 
formulation of the collision rules its physical interpretation as a degree of freedom is now more flexible. For example, 
one may think of k as being related to some kind of lattice modes. 

To ensure that the system is time-reversible, we let the forward baker act if < 7 < t^/^, and its inverse if 
— 7r/2 < 7 < 0. The angle 7' always goes to the respective other side of the normal, that means 7' has the opposite 
sign of 7, as shown in Fig.y. To avoid a symmetry breaking in a possible nonequilibrium situation we alternate the 
assignment of the forward and backward baker, that is, if int(/9 • 10*) is even (odd) we take the forward (backward) 
baker for 7 > and vice versa the corresponding backward (forward) baker for 7 < 0. Ideally, the alternation should 
be done in infinitely fine steps, which is not feasible in computer simulations. 

2. Relation between map density and time continuous density 

By investigating the dynamics of the collision process through a baker map, one is actually considering the Poincare 
section of the velocity of the particle at the moment of the collision. We denote the corresponding probability density 
of the moving particle as the map density Qmap{v). Qmap{v) can be written in discretized form as Qmap{vi) = 



*Note that the Poincare or Birkhoff mapping between two collisions is just area preserving in the Birkhoff coordinates 
(/3,sin(7)). 
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(. . . , Cvjco, ■ ■ .), where c^- is the number of coUisions after which the particle has the velocity Vi and Cq is the total 
number of collisions. 

One may establish a relation between Qmap{v) and the time continuous probability density q{v)^ where v is measured 
at any time interval dt. This is the relevant quantity to check for a microcanonical distribution. Notice that we can 
write in the same way as before Q{vi) = (..., i^./t, . . .), where is the total time during which the particle has the 
velocity w^, and t is the total time. 

We introduce now the mean time of flight between two collisions < i > by 

< <>- - (2) 

Co 

and the mean time of flight between two collisions < t >„. when the particle has the velocity Vi by 

<<>.,= ^ = ^^. (3) 

< s > is the collision length, which is expected to be independent of v. To compute a value of givi), Eq. (H) and 
Eq. can be combined to 

t < t > Vi Cq 

Eq. is valid for all i implying that we get the following relation between the time continuous density g(v) and the 
map density g„iap{v), 

^ 9maM < S > ^ ^^^^^ emaAv) (5) 
V < t > V 

where the constant is determined by normalization. By assuming that the coupling of the collision rules of Eq. (Q) to 
the specific geometry of the Lorentz gas does not yield an invariant map density being different from the one of the 
baker map, we can now calculate g{v) for our model: Inserting Qbaker = 1 in Eq. (^) results in g{v) — const. /v. This 
expression is not normalizable and does not correspond to the correct result related to a microcanonical density. The 
way out of this difficulty is presented in the next section. 

3. Getting a microcanonical density 

Having shown in the last subsection that the collision rules as described by the baker map are not sufficient to 
get a microcanonical probability density, we now amend the definition of the collision rules. We do this by including 
an additional transformation Y which is constructed in the following way: Our system has three degrees of freedom 
{vx, Vy,k), and the total energy E is conserved, 2E — v^'^ +Vy'^ + k'^. The dynamics is microcanonical if the probability 
density is equidistributed on a three-dimensional sphere, 

g3{vx,Vy,k) = —^S{2E - Vx'^ - Vy'^ - k^) . (6) 
In the Appendix we calculate the reduced density for one or two degrees of freedom of a d-dimensional microcanonical 



system. For d = 3, Eq. (All) and Eq. (A17) lead to 



With Eq. (|) the map density corresponding to Eq. reads 

To get the reduced microcanonical probability densities given by Eqs. (Q) and (^) for the velocity of the particle in 
our model we thus have to redefine the baker variable yb- Conservation of probability 
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Qbaker{yb)dyb = Qmap{v)dv 



(10) 



yields 

Vb = Y3,{v) = - — \J2E - w2 + - arcsin— L= (11) 

with < u < \J2E , < ¥3(1;) < 1. The inverse transformation v — Y^^^i^yb) exists because Yj, is monotonous. With 
Xh = ^(7) = sin I7I, we can summarize the cohision rules to 

(7', v') = {X-\Y^^) o M o (X(7), Y^{v)) . (12) 

fc' being obtained from energy conservation, k' — \/ 2E — v'^. Fig. ^ shows the probability densities for the three- 
dimensional system resulting from numerical simulations at i? = 0.5. g3{(3) = gs^a) = l/27r are uniform as expected. 
^3(^)7 Qai^x) and Qaivy) are in exact agreement with the reduced microcanonical densities, that is, gs^v) corresponds 
to Eq. (^ and gsivx), g^ivy) correspond to Eq. (0). Moreover, numerical simulations show that the trajectory of the 
particle covers the Lorentz gas cell uniformly in configuration space. Beside the equidistribution of gsiP), this is a 
further check of the ergodic behavior of our system. 



C. Arbitrarily many degrees of freedom on the disk 

We have considered an energy transfer between particle and disk for the case when the disk is equipped with one 
degree of freedom. We now further modify the dynamics by pretending that the disk has arbitrarily many degrees 
of freedom k = {ki, . . . ,kd-2), entailing that in total the system has d degrees of freedom. We do not deal with 
the individual components of fc, because the detailed dynamics of k is not relevant for our purpose. Instead, we 
consider only the absolute value |fc|. The description of our model as a dynamical system is therefore still based on 
having only three relevant variables for the velocities, whereas the corresponding statistical physical situation involves 
a microcanonical probability density of a d-dimensional system, and thus implicitly mimics the situation of having d 
physical degrees of freedom. 

The consideration of such additional degrees of freedom requires a modification of the collision rules to get correctly 
the microcanonical probability density corresponding to this d-dimensional system. In particular, we have to redefine 
the transformation Y, which can be done by the same method as for the three dimensional system above: First, 
we c alcul ate the reduced densities gd{vx) and gd{v) of the d-dimensional energy hypersphere. They are given by 
Eq. (|All| ) and Eq. (|A17| ) in the Appendix, 

gM = - (13) 

VttI (— ) {2E) 2 

and 

gd{v) = -^—^v[2E ~ v^)^ . (14) 
{2E)^ 

Inserting Eq. (^4|) into Eq. (||) leads to the corresponding map density, and with Eq. ( p^ we can calculate Yd- For 
arbitrary even d it is given by 

I 



Yd = \~v{2E-v')— + y \ 2 \^ '-v''+'{2E) — -'\. (15) 

r(|)r(^)(2i?)^(d-2) ^ t^\^)2^ + l ^ ^ I 

For all d this expression is monotonous, and thus its inverse always exists. For odd d> ZYd can be calculated in the 
same way, but cannot be written down in a closed form. Inserting Eq. ( p^ ) in Eq. (|l2|) yields the full collision rules of 
our model for arbitrary even d. \k'\ is obtained from energy conservation. Fig. 3| shows computer simulation results 
of Pd(v) and gd{vx) for d = 3, 4 and 6. They are in exact agreement with Eqs. (13|) and (|lj). 

Eq. (13) was already written by Maxwell and Boltzmann ||7|j7^, however, for their calculation they took a different 



starting point: a set of n particles, any particle with two degrees of freedom, is moving according to Hamilton's 
equations of motion. Using properties of the canonical transformation they derived the reduced microcanonical 
probability densities Eqs. ( p^ ) and (p^). Notice that starting the calculation of the reduced densities according to our 
reasoning, but including also momentum conservation, reduces the effort of the calculation drastically, see Ref . pS] . 
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D. The disk as a thermal reservoir 



We now c onsid er the hmit d — > oo. Using equipartitioning of energy E — dT/2 with — 1, Eq. (All) reduces 



then to Eq. ( |A14| ) 



1 ^ 

Q^{vx) = -^=e-5t , (16) 

which is the canonical distribution for the moving particle. The disk now acts as a thermal reservoir, yet, the whole 
system remains microcanonical. In the same way the limiting form of Eq. (A17), see Eq. ( Al^ ), reads 



1 

g^{v) = j;ve . (17) 



As for the transformation Yoo, it is given by 



Y^{v)^-^—ve^+eTf{^). (18) 



As before, Yoo is monotonous, and thus the inverse exists. Notice that in all these expressions the temperature 
T appears instead of the previous total energy of the system and serves as a free parameter. 

Fig. ^ shows the probability densities for d ^ oo resulting from computer simulations: QooiP) and £>oo(a) are 
uniform, as expected. goo{v) is in exact agreement with Eq. ( |l7|) and goo{vx) and gooivy) are in exact agreement with 
Eq. ®. 

We note that combining Eq. (11^) with the respective equilibrium distribution for the angle of incidence goo (7) — cos 7 
and applying Eq. (^ yields a map density gmapiliv) in local polar coordinates which is identical to the stochastic 
boundary conditions as given, e.g., in Ref. Js^ . This relation of our method to stochastic boundary conditions is 
discussed in more detail in Ref. | |7^ , where it has been used as an alternative starting point to define the deterministic 
and time-reversible thermostat constructed above. 

III. NONEQUILIBRIUM, STEADY STATES, AND TRANSPORT UNDER AN EXTERNAL ELECTRIC 

FIELD 

In the previous section we have extended the dynamics of a periodic Lorentz gas in equilibrium by defining new 
collision rules for the particle allowing for an energy transfer between particle and disk at each collision such that the 
dynamics is deterministic, time reversible, and yields a microcanonical probability density. In this section we analyze 
the behavior of the model under nonequilibrium conditions, associated with the presence of an external electric field 
£ applied parallel to the x-axis. Taking the collision rules as defined in equilibrium, we study the structure of g(v) in 
nonequilibrium and the associated transport properties. 

Fig. ^ shows the time evolution of twice the average kinetic energy, < >, for an ensemble of moving particles 
and for different d's. As can be seen, the particle energy grows continuously with time for finite d, but fluctuates 
around a constant mean value as d — > 00. This can be understood as follows. In equilibrium the energy transfer 
ensures equipartitioning of energy between all degrees of freedom. In the presence of a field the energy of the particle 
grows during the free flight, and as a consequence the particle has on the average a surplus of energy at a collision 
in comparison to the disk. The energy transfer then counteracts this surplus of energy of the particle, because 
equipartitioning of energy is still built into the collision rules. Since there is no other source of dissipation in our 
model than the transfer of energy onto the disk, the energy of the particle must eventually grow for finite dimension, 
while the growth rate decreases by increasing d, because by increasing d more energy can be stored onto the disk. In 
the limit of d — > 00 we obtain a constant average kinetic energy, since then the disk acts as a thermal reservoir with 
infinitely many degrees of freedom, which means that our system is thermostated. Still, the system is time reversible. 



This interesting fact is reflected in the special functional form of the transformation Yoc{v) in Eq. (18). 

In the following we investigate the nonequilibrium steady state of the model for c? — > 00 in more detail. The mean 
kinetic energy in the comoving frame < > — < w^; >^ of the moving particle in comparison to the one obtained from 
the equipartition theorem, < v'^ > — < Vx 2T, is presented in Fig. |6[ For small e, or for large temperatures, the 
curves approach the equipartitioning values. However, for general s and T there is a systematic difference between 
the measured temperature in the simulations and the parametric temperature as it appears in the collision rules. 
The reason for this difference can be found in our approach to define the collision rules in an equilibrium situation: 
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First, Eq. (|5|), which is one step in the derivation of the transformation ^00, is based on having a constant particle 
vefocity between two colhsions, which is not correct anymore by applying an electric field. And second, we always 
try to transform the system onto an equilibrium distribution by assuming that this will still be a good approximation 
for a system slightly perturbed out of equilibrium. We note that we do not have any better choice, since we do not 
know the correct noncquilibrium distribution for the Lorentz gas. As a consequence, the parametric temperature of 
the disk, which refers to an equilibrium distribution, and the measured kinetic temperature of the moving particle 
in noncquilibrium do not agree. One way out of this problem would be to redefine a proper kinetic temperature for 
the reservoir. This can in fact be done by employing the average value of the variances of the in- and outgoing map 
densities, as is discussed in detail in Ref. [Q. With respect to such a new temperature definition for the reservoir, we 
would expect to have equipartitioning of energy at least for sufficiently small field strengths when the temperature is 
sufficiently high. Fig. ^ shows that for low T the variance approaches finite values not equal to zero for non- vanishing 
field strengths. The reason for this behavior is that between collisions the electric field accelerates the particle, 
whereas the thermostat acts only at the collisions. Thus, we may expect that in this region there will always remain 
a discrepancy between the measured temperatures of particle and reservoir. 

Fig. I?! depicts the probability densities for d ^ 00 and e = 0.5. The external field leads to strong deviations from 
the equilibrium probability densities: On a very coarse scale, Qoo{P) shows the existence of a global maximum opposite 
to the field direction at /3 = tt. However, there appear six strong pronounced local maxima which are related to the 
specific geometry of our model. Correspondingly, £Ico(q;) shows a global maximum parallel to the field direction at 
a = 0, but there also exist four local maxima on a finer scale. Qooi'v)^ Qooivx) and Qoa{vy) are also clearly modified 
by the field, while remaining close to the functional form of the equilibrium probability densities. Especially, Qoo{vx) 
is shifted to positive Vj^-values. This indicates the existence of a current parallel to the field direction. 

The changes in the probability densities by increasing the field strength in reference to the equilibrium solutions 
are presented in Fig. ^. The mean value of (?oo(t^a;) grows with e, whereas Qoa{vy) remains symmetric around 0. 

The conductivity a =< > /e, shown in Fig. |9|, is a globally decreasing function of e. This demonstrates that 
for the field strengths considered in the figure we are already in a highly nonlinear regime. We furthermore note that 
there exist irregularities in ct on a finer scale, which are beyond our numerical error estimates and indicate a deviation 
of a from a simple functional dependence on e. This curve should qualitatively be compared to the conductivity as 
obtained for the Gaussian thermostated periodic Lorentz gas[] ||l9|,^,^ : Fig. 2 in Ref. shows a globally decreasing 
conductivity on a coarse scale as well, however, its irregularities on a fine scale are much more pronounced and clearly 
nonmonotonous in e. Whether the irregularities in our Fig. ^are in fact monotonous or nonmonotonous on a finer scale 
cannot be decided on the basis of our numerical data. Unfortunately, it is not clear how to compare the conductivities 
of these two models quantitatively, since the choice of temperature in the Gaussian model is somewhat ambiguous by 
a factor of two [|2|,||. 

To understand why no linear response is seen in Fig. |^, we attempt an estimate of the expected range of validity 
of the linear response regime by applying the simple heuristic argument suggested, e.g., in Refs. and further 

references therein. The reasoning given in these references may actually be considered as the standard response to 
the famous van Kampen objection against linear response ||75| , since it amounts to properly modifying his original 
argument. For the periodic Lorentz gas it is stipulated that for having linear response it is sufficient for the field 
strength to fulfill e ^ '^/'''coii' where TcoU is the average mean- free time between collisions. We have computed this 
bound for the Gaussian thermostated Lorentz gas as well as for our model. For the Gaussian thermostat, TcoIi has 
already been obtained from computer simulations at different densities of scatterers in Ref. |^^. For the density 
corresponding to the gap size w — 0.2361 the upper bound for linear response is then approximately e <C 5.6. 
However, numerical results for the conductivity indicate no linear response down to e ~ 0.05 ||67|| . For a higher 
density corresponding to w = 0.129, the conductivity in Fig. 9 of Ref. and the respective data values may be 
compatible with a linear response-like behavior at most up to a field strength of approximately e ~ 0.9. Here, an 
upper bound for linear response yields e <C 13. Finally, for w = 0.0076 the respective upper bound is e ^ 70. 
The corresponding conductivity appears to be well-behaved up to s ~ 2.5, although for this very high density the 
numerical values are not very precise. For our model, we calculated as an upper bound s <C 4.45, however, down to 
a field strength of e = 0.05 we cannot detect any linear response. 

We conclude that in the light of the numerical simulations the bound for linear response discussed above does 
not appear to serve as a very reliable reference for the driven periodic Lorentz gas. Other restrictions, which may 
make this bound more precise, have been discussed in Ref. [^3|| , without leading, however, to an improved quantitative 
formulation. Thus, although the existence of a linear response has been proven for the Gaussian thermostated periodic 



tNote that in Refs. PBl the angle between the field direction and the x-achses has been chosen equal to ^. 
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Lorentz gas with an external field in Refs. |^,^, having a reasonable estimate for the range of validity of linear 
response in this system remains, to our knowledge, an open question. In particular, we note that up to now for 
smaller densities like w = 0.2361 such a regime has not been detected in computer simulations. We furthermore note 
that analogous problems, namely the existence of irregularities on a fine scale in parameter-dependent deterministic 
transport coefficients as well as apparently the complete lack of a linear response regime have already been encountered 
in more simple low-dimensional models of chaotic transport. These are so-called multibaker maps, which are believed 
to represent certain essential dynamical features of Lorentz gas models [ |56| , [76| , [77t , and associated one-dimensional 
maps. The currents in such systems have been computed by numerically exact methods with respect to varying a 
bias, and they turned out to be fractal functions in the bias 78-8^], implying also the non-existence of a regime of 



linear response down to extremely small field strengths. Such properties may be related to the low dimensionality of 
these systems, which is shared by the periodic Lorentz gas. 

Fig. (a) shows the Poincare section of (sin(7), (3) at the moment of the collision obtained for our model with the 
baker map. It displays a highly variable phase space density exhibiting a structure which is qualitatively analogous 
to the one of the fractal attractor found for the Gaussian thermostated Lorentz gas ]l9| , p^ , p6l|27y43| ] . The existence 
of such an attractor for our model is due to the fact that the phase space variables in Fig. |10| (a) only reflect the 
dissipative dynamics of the moving particle, whereas the corresponding complementary dynamics of the reservoir is 
completely left out. However, the dynamics of the single moving particle, which represents here a subsystem projected 
out of the full system (particle plus reservoir), is certainly not phase space volume preserving, because there is an 
average energy transfer to the reservoir. This does not contradict the existence of volume preservation, and thus a 
related Hamiltonian character, of the full system (particle plus reservoir) if one takes the complete dynamics of the full 
system appropriately into account. This reasoning can be made more precise by computing the Jacobian determinant 
for the full system (particle plus reservoir) at a collision by systematically varying the number of degrees of freedom 
associated to the reservoir [ |8l|] . 

Figure |l^ (b) shows the Poincare section for a modification of our model where we have replaced the baker map 
by a random number generator, thus modeling in a way stochastic boundary conditions. This stochastic model 
apparently leads to a smooth nonuniform density. The existence of smooth versus singular measures in thermostated 
dynamical systems like the driven Lorentz gas has been vividly discussed in the recent literature |l^,^,^,^: As 
pointed out in the introduction of this paper, the existence of singular nonequilibrium measures is an essential feature 
of Gaussian and Nose-Hoover thermostated systems. On the other hand, it has been proven that for a specific 
system in a nonequilibrium situation created by stochastic boundaries the corresponding measure is smooth [p^ , 
and it has been argued that this should be considered as a general feature of systems which are thermostated by 
stochastic boundaries |^^. Our results of Figs. |l^ (a) and (b) suggest that for the periodic Lorentz gas there is a 
clear distinction between deterministic boundaries producing a singular measure, and stochastic boundaries creating a 
smooth measure. We note that in Ref. a combination of deterministic and stochastic boundaries has been applied 
to the driven periodic Lorentz gas, leading to the observation that the fractal attractor is apparently very stable 



with respect to the random perturbations induced by the stochastic part of the boundaries. In Ref. |82| computer 
simulations appear to indicate the existence of a slightly singular measure for a system with stochastic boundaries. 
How the nonequilibrium stationary measure looks like for a general many-particle systems with stochastic boundaries 
remains thus an open question up to now. Our approach of thermostating may help to bridge the gap, because we can 
alternatively produce deterministic or stochastic boundaries simply by replacing the reversible deterministic map by 
a random number generator, thereby either preserving or completely destroying any dynamical correlations induced 
by the boundaries. 

Fig. |l^ (c) shows the Poincare plot for /3 — f{e) a.t T — 0.5. For small e we observe a uniform phase space density 
which is in agreement with the histogram £ioo(/3) in Fig. |^. By increasing e the density shows several maxima but 
remains phase space filling. This behavior refiects the dynamics depicted in the histogram Poo(/3) of Fig. 0. We 
could neither detect a contraction of the attractor onto periodic orbits, nor a breakdown of ergodicity for higher field 
strengths, or eventually the existence of a so-called creeping orbit for very high field strengths^. Similar results have 
been obtained for other choices of Poincare sections in phase space. This is in contrast to what has been found for 
the Gaussian isokinetic thermostat |^Jl^,^,^,^, suggesting that the complicated scenario found in the Gaussian 
model is rather a property induced specifically by the Gaussian thermostating mechanism. 



■'■If the alternation of forward and backward baker as introduced in section 
windows, which disappear in the limit of infinitely fine steps. 



[I B 1 is not fine enough one observes periodic 
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IV. CONCLUSIONS 



We have presented a new deterministic thermostating mechanism for the periodic Lorentz gas. Our mechanism is 
based on modeUng the energy transfer between the particle and the disk at each coUision instead of using a momentum- 
dependent friction coefficient. We have defined the coUision rules for the particle such that they are time reversible, 
deterministic, and yield a microcanonical probability density in equilibrium. We have shown how our mechanism 
could be applied for a disk equipped with arbitrarily many degrees of freedom. We investigated the behavior of the 
model under nonequlibrium conditions by applying an external electric field. In the limit d — > oo, where the disk 
acts as a thermal reservoir, we obtained a nonequilibrium steady state with a constant average kinetic energy of the 
particle implying that we have successfully thermostated our system. 

With respect to the parametric temperature contained in our collision rules the dynamics of our model nearly 
fulfills equipartitioning of energy for small field strengths, however, by increasing the field strength the deviation from 
equipartitioning grows. This appears to be a problem of a definition of temperature in nonequilibrium and could 
probably be avoided by using an appropriately defined temperature for the disk. 

The results obtained from our model in nonequilibrium have been compared to the results of the Gaussian ther- 
mostated periodic Lorentz gas. It turned out that the attractor associated to the model is qualitatively similar to 
the fractal attractor of the Gaussian thermostated periodic Lorentz gas. However, we could not find any breakdown 
of ergodicity or a coUaps of the attractor onto periodic orbits for higher field strengths in our model. In both mod- 
els the conductivity is a nonlinear decreasing function with increasing field strength with more or less pronounced 
irregularities on a finer scale. 

Our thermostating mechanism yields the canonical probability density for the moving particle in equilibrium and 
keeps the kinetic energy of the particle constant on average in nonequilibrium. This is similar to what is achieved 
by applying a Nose-Hoover thermostat. However, to our knowledge a Nose-Hoover thermostated Lorentz gas has 
not yet been studied so far. Work along these lines and comparison with the results of our model is currently under 
investigation fs^ . 

We note that the thermostating mechanism presented in this work has recently been successfully applied to simulate 
a heat and shear flow of an interacting many-particle system of hard disks |7j]. As a main result, it has been found 
that in general there exists no identity between phase space contraction and entropy production in this system if it is 
thermostated by our method. We would expect the same result to hold for the thermostated Lorentz gas as described 
in this paper. It would also be interesting to investigate the validity of fluctuation theorems p^ , |8^ for our Lorentz 
gas model. 

Further work along the lines of this paper should also aim at calculating the Lyapunov exponents, using for 
instance the method of Dellago and Posch ji^js^,^ . This would enable to check for the validity of the expression 
of the conductivity in terms of the sum of Lyapunov exponents, as obtained for conventionally thermostated systems 
p^ , ^ 32 Moreover, it could be investigated whether there exists a conjugate pairing rule for the Lyapunov 

exponents in our model ||7|,^o|]. In particular, the knowledge of the Lyapunov exponents would allow an estimation of 
the fractal dimension of the attractor. 

Finally, collision processes in granular media might be advantageously modeled by applying our formalism of energy 
transfer. This could be done by associating finitely many degrees of freedom to a scatterer according to our method, 
thus rendering the collision process inelastic. This would provide an interesting alternative to describing collision 
processes by velocity dependent restitution coefficients, as it has been done previously 
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APPENDIX A: CALCULATION OF THE REDUCED DENSITIES 



In this section we calculate the reduced, or projected, densities for one and two degrees of freedom of a d-dimensional 
microcanonical system. We use the following notation: general momentum space coordinates are denoted by xi, . . . , 
or by X, and Wi is the absolute value of a vector with i , i €i N, degrees of freedom. 

The probability density of a d-dimensional microcanonical system is an equidistribution on a d-dimensional hyper- 
sphere, 

Qdixi, ...,x,)^ ^M^Si2E ^xl-...-xl). (Al) 
With the definition of rf-dimensional spherical coordinates 
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xi = V2E simpd-1 sin'ipd-2 ■ ■ ■ sin V'2 sinipi 
X2 — V2E sin V-'d-i sin 4'd~2 ■ ■ ■ sin 'ip2 cos ipi 
X3 — V2E sin V'd- 1 sin ■(/'d- 2 ■ • • cos V'2 



Xd-l 
Xd 



^2£'sin T^d-i cos i>d-2 
^2Ecos'ipd-i , 



(A2) 



where -01 G [0, Stt) and tjjn € [0, tt), n > 2, we can transform Eq. ( Al) onto these coordinates according to 

QdiV^E, "01, • ■ • , ^d-i) = Jd{xi, . . . V2E, V'l, ■ • ■ , '0<i-i)£'d(a;i, . . . , (A3) 

where Jd = \{dxi . . . dxd) / {d^/2Edipd-i ■ ■ ■ dtpi)\ is the Jacobian determinant. For Jd, d > 3, we obtain the following 
recursion relation, 



Jd = 



sinipd-ian V2E coatpd-ian sin'ijjd-iai2 
sinV'd-ia2i v^cosV'd-ia2i sin'0d_iQ;22 



sin tpd-iai{d~-i) 
sinV'rf_iQ!2(<i-i) 



sinV'd-ia(d-i)i V2Ecos^pd-ia(^d-i)i sinipd~ia{d~i)2 •■• sin 

cosi/'d-i — V2^sin-0d-i ... 



(A4) 



where Ujk are the matrix elements of the Jacobian Jd-i for d — I dimensions. An expansion after the d-th row yields 



Jd = (cos?/'<i-i|f^i| + V2Esm-)pd-i\^2\ 



(A5) 



Any column in the matrices f2i and fl2 is equivalent to the corresponding column in Jd-i multiplied by a factor. In 
the first column of Oi (derivative after ipd-i) the factor is V2£'cos ipd-i, in the other columns of Qi the factor is 
s'mipd-i- In all columns of $^2 the factor is sinipd-i- Factoring out leads to 



Jd = V2Esm'^-^'tljd-iJd^i 
Proceeding for the Jacobians Jd-i ■ ■ ■ J2 in the same way gives 

d-i 



Jd = {2E)—l[sm"-'ij^ 



(A6) 



(A7) 



n=2 



Inserting Eq. (A7) into Eq. ( [A3| ) and taking into account that the radius appearing in the spherical coordinates is 
a constant equal to V2E, according to 6{2E — — ... — x^) = 1 leads to the microcanonical probability density in 
spherical coordinates 



rf-i 



z(V'i,...,^d-i)=^j^n«in""'^" 

27r2 



(A8) 



In Eq. (A8) the angles ipn appear in product form. This allows to infer from this equation directly the probability 
density for the angle 'ipn. 



Normalizing g{ipn), n >2, leads to 



e(V„) =iVsin"-iVr, 



(A9) 



(AlO) 



and g{ipi) — l/27r. Now we can calculate the probability density for a variable representing exactly one degree of 
freedom by using the last hne of Eq. (A2), 
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Qdixd) = g{ipd~i)\ 



gd[xd) 



dxd 

r(|) 1 



— {2E~xl) 



(All) 



Eq. (All) is valid for all 1 < n < d. To simplify th e no tation we replace the variable Xd by x in the following. 
Inserting equipartitioning of energy E = Td/2 into Eq. (All) and performing the limit d oo leads to 



lim (1 -) 2 = e 2T 

d^oc Td 



for the last term and, by applying Stirlings formula, to 

r(f 



for the prefactor. Finally, we get 



lim — 3 — 



Qooix) = lim gd{x) = — f 
d^co v27rr 



(A12) 



(A13) 



(A14) 



The probability density gdiwd-i) with Wd-i 
lated straightforward to 



Qd{wd-i) = 2g{-)pd-i)\ 
Qd{wd~i) 



dtp. 



d-l 



dwd-1 

r(i) 2 



.T^ j^ — \JlE — = v2^ sin V'd-i can be calcu- 

(A15) 



d-2 



^/^^(^) (2i?)^ ^ 



Note that the angles ipd^i^^ < V"!;-! < 7 r/2 a nd 7r/2 < tt — i/irf-i < tt are mapped onto the same value w^-i, Wd-\ > 0, 
which leads to the factor of two in Eq. ( A15 ). 

We proceed by considering the variable Wd-2 — \Jx\ + x\ + . . . + x'^^_^ — Wd-i svLV%lid-2- Together with the variable 

Xd-i = Wd-i costJjd-2 of Eq. ( A2) a two-dimensional set of variables {wd-2, Xd-^i) can be defined, where (w^-i, '0d-2) 
are the corresponding polar coordinates. On these grounds the probability density Qd{wd-2) can be calculated to 



Qd{Xd-\,Wd^2) = Qd{wd-l)Q{ll;d~2)\ 

d-2 



d{wd-i,ipd-2) I 

d{xd^l,Wd^2) 
^d-2 



n{2E)- 



2^ - ^i-2 ~ ^l-i 



Integration over Xd-i yields 



Qd{wd-2) 



d-2 



{2E) 



,,<i-3 

d-2 ^d-2 



(A16) 



With W2 = a/2-E' — ""^^-2 the probability density for two degrees of freedom Qd{w2), 



Qd{w2) 



d-2 



{2E) 



—W2{2E-wi)- 



Inserting E — Td/2 and performing the limit d — > oo yields 

£-00(^2) = — e 



(A17) 



(A18) 



11 



[1] W. G. Hoover and W. T. Ashurst, in Theoretical Chemistry, edited by H. Eyring and D. Henderson (Academic, New York, 
1975), Vol. 1. 

J. M. Haile and S. Gupta, J. Chem. Phys. 79, 3067 (1983). 

M. D. Allen and D. J. Tildesley, Computer simulation of liquids (Clarendon Press, Oxford, 1987). 
D. J. Evans and G. P. Morriss, Statistical Mechanics of Nonequilibrium Liquids (Academic Press, London, 1990). 
W. G. Hoover, Computational statistical mechanics (Elsevier, Amsterdam, 1991). 

S. Hess, in Computational physics, edited by K. Hoffmann and M. Schreiber (Springer, Berlin, 1996), pp. 268-293. 
G. P. Morriss and C. P. Dettmann, Chaos 8, 321 (1998). 

W. G. Hoover, A. J. C. Ladd, and B. Moran, Phys. Rev. Lett. 48, 1818 (1982). 
D. J. Evans, J. Chem. Phys. 78, 3297 (1983). 
D. J. Evans et al, Phys. Rev. A 28, 1016 (1983). 
S. Nose, Mol. Phys. 52, 255 (1984). 
S. Nose, J. Chem. Phys. 81, 511 (1984). 
W. G. Hoover, Phys. Rev. A 31, 1695 (1985). 
J. L. Lebowitz and H. Spohn, J. Stat. Phys. 19, 633 (1978). 
A. Tenenbaum, G. Ciccotti, and R. Gallico, Phys. Rev. A 25, 2778 (1982). 
S. Goldstein, C. Kipnis, and N. laniro, J. Stat. Phy s. 41, 915 (1985). 



R. Klages, K. Rateitschak, and G. Nicolis, preprint |shao-dyn/981202l| (1998). 
D. J. Evans and B. L. Holian, Phys. Rev. A 83, 4069 (1985). 
B. Moran and W. G. Hoover, J. Stat. Phys. 48, 709 (1987). 

B. L. Hohan, W. G. Hoover, and H. A. Posch, Phys. Rev. Lett. 59, 10 (1987). 
W. G. Hoover, Phys. Rev. A 37, 252 (1988). 

G. P. Morriss, Phys. Lett. A 122, 236 (1987). 

G. P. Morriss, Phys. Lett. A 134, 307 (1989). 

W. G. Hoover and B. Moran, Phys. Rev. A 40, 5319 (1989). 

G. P. Morriss, Phys. Rev. A 39, 4811 (1989). 

W. G. Hoover and H. A. Posch, Chaos 8, 366 (1998). 

W. G. Hoover and H. A. Posch, Phys. Lett. A 246, 247 (1998). 

G. Eyink and J. Lebowitz, in Microscopic simulations of complex hydrodynamic phenomena, Vol. 292 of NATO ASI Series 
B: Physics, edited by M. Mareschal and B. L. Holian (Plenum Press, New York, 1992), pp. 323-326. 

H. A. Posch and W. G. Hoover, Phys. Lett. A 123, 227 (1987). 
H. A. Posch and W. G. Hoover, Phys. Rev. A 38, 473 (1988). 
H. A. Posch and W. G. Hoover, Phys. Rev. A 39, 2175 (1989). 

N. L. Chernov, C. L. Eyink, J. L. Lebowitz, and Y. G. Sinai, Phys. Rev. Lett. 70, 2209 (1993). 

N. L. Chernov, C. L. Eyink, J. L. Lebowitz, and Y. G. Sinai, Comm. Math. Phys. 154, 569 (1993). 

N. L Chernov and J. L. Lebowitz, Phys. Rev. Lett. 75, 2831 (1995). 

T. Tel, J. VoUmer, and W. Breymann, Europhys. Lett. 35, 659 (1996). 

J. VoUmer, T. Tel, and W. Breymann, Phys. Rev. Lett. 79, 2759 (1997). 

N. L Chernov and J. L. Lebowitz, J. Stat. Phys. 86, 953 (1997). 

G. Gallavotti and E. G. D. Cohen, Phys. Rev. Lett. 74, 2694 (1995). 

D. Ruelle, J. Stat. Phys. 85, 1 (1996). 

D. J. Evans, E. G. D. Cohen, and G. P. Morris, Phys. Rev. A 42, 5990 (1990). 
W. N. Vance, Phys. Rev. Lett. 69, 1356 (1992). 

A. Baranyai, D. J. Evans, and E. G. D. Cohen, J. Stat. Phys. 70, 1085 (1993). 

C. Dellago, L. Glatz, and H. A. Posch, Phys. Rev. E 52, 4817 (1995). 

J. R. Dorfman, An introduction to chaos in nonequilibrium statistical mechanics (Cambridge University Press, Cambrigde, 
1999). 

Chaos and Irreversibility, Vol. 8 of Chaos, edited by T. Tel, P. Gaspard, and G. Nicolis (American Institute of Physics, 
College Park, 1998). 

Microscopic simulations of complex hydrodynamic phenomena. Vol. 292 of NATO ASI Series B: Physics, edited by M. 
Mareschal and B. L. Holian (Plenum Press, New York, 1992). 

The microscopic approach to complexity in non- equilibrium molecular simulations. Vol. 240 of Physica A, edited by M. 

Mareschal (Elsevier, Amsterdam, 1997). 

W. G. Hoover et al, Phys. Lett. A 133, 114 (1988). 

C. P. Dettmann and G. P. Morriss, Phys. Rev. E 54, 2495 (1996). 

C. P. Dettmann and G. P. Morriss, Phys. Rev. E 55, 3693 (1997). 

P. Choquard, Chaos 8, 350 (1998). 

C. Dellago and H. A. Posch, J. Stat. Phys. 88, 825 (1997). 
P. Gaspard and G. Nicolis, Phys. Rev. Lett. 65, 1693 (1990). 
P. Gaspard and J. R. Dorfman, Phys. Rev. E 52, 3525 (1995). 



12 



[55] J. R. Dorfman and P. Gaspard, Phys. Rev. E 51, 28 (1995). 

[56] P. Gaspard, Chaos, Scattering, and Statistical Mechanics (Cambridge University Press, Cambridge, 1998). 
[57] H. A. Lorentz, Proc. Amst. Acad. 438 (1905). 

[58] L. A. Bunimovich and Y. G. Sinai, Commun. Math. Phys. 78, 479 (1981). 
[59] P. Cvitanovic, P. Gaspard, and T. Schreiber, Chaos 2, 85 (1992). 
[60] J. Machta and R. Zwanzig, Phys. Rev. Lett. 50, 1959 (1983). 
[61] P. Gaspard, Chaos 3, 427 (1993). 

[62] P. Gaspard and F. Baras, Phys. Rev. E 51, 5333 (1994). 

[63] G. P. Morriss and L. Rondoni, J. Stat. Phys. 75, 553 (1994). 

[64] P. Gaspard, Phys. Rev. E 53, 4379 (1996). 

[65] H. Matsuoka and R. Martin, J. Stat. Phys 88, 81 (1997). 

[66] J. Lloyd, L. Rondoni, and G. P. Morriss, Phys. Rev. E 50, 3416 (1994). 

[67] J. Lloyd, M. Niemeyer, L. Rondoni, and G. P. Morriss, Chaos 5, 536 (1995). 

[68] C. P. Dettmann and G. P. Morriss, Phys. Rev. E 54, 4782 (1996). 

[69] C. P. Dettmann and G. P. Morriss, Phys. Rev. Lett. 78, 4201 (1997). 

[70] K. Rateitschak and R. Klages (unpublished). 

[71] J. C. Maxwell, Cam. Phil. Trans. 12, 547 (1879). 

[72] L. Boltzmann, in Wissenschaftliche Abhandlungen von L. Boltzmann, edited by F. Hasenohrl (J. A. Barth Verlag, Leipzig, 

1909), Vol. 2, Chap. 63. 
[73] L. J. Milanovic, H. A. Posch, and W. Thirring, Phys. Rev. E 57, 2763 (1998). 
[74] C. Wagner, R. Klages, and G. Nicolis, Phys. Rev. E 60, (in press) (1999). 
[75] N. van Kampen, Physica Norvegica 5, 279 (1971). 
[76] P. Gaspard and R. Klages, Chaos 8, 409 (1998). 
[77] W. Breymann, T. Tel, and J. Vollmer, Chaos 8, 396 (1998). 
[78] J. Groeneveld, priv. commun. 

[79] R. Klages and J. Groeneveld, Verhandl. DPG VI, 646 (1998). 
[80] R. Klages, Verhandl. DPG VI, 678 (1999). 
[81] R. Klages (unpublished). 

[82] H. A. Posch and W. G. Hoover, Phys. Rev. E 58, 4344 (1998). 

[83] K. Rateitschak, R. Klages, W. Hoo ver, and G. Nicolis ( unpublished). 

[84] J. Lebowitz and H. Spohn, preprint |cond-mat/9811220 . 

[85] C. Dellago and H. A. Posch, Phys. Rev. E 52, 2401 (1995). 

[86] C. Dellago, H. A. Posch, and W. G. Hoover, Phys. Rev. E 53, 1485 (1996). 

[87] N. Brillantov, F. Spahn, J.-M. Hertzsch, and T. Poschel, Phys. Rev. E 53, 5382 (1996). 

[88] T. Poschel and T. Schwager, Phys. Rev. Lett. 80, 5708 (1997). 

[89] T. Schwager and T. Poschel, Phys. Rev. E 57, 650 (1998). 

[90] G. Ciese and A. Zippelius, Phys. Rev. E 54, 4828 (1996). 

[91] T. Aspelmeier, G. Giese, and A. Zippelius, Phys. Rev. E 57, 857 (1998). 



13 




14 




FIG. 3. Probability densities for d = 3 (solid curve), 4 (dashed curve) and 6 (long dashed curve) at E = 0.5 




FIG. 4. Probability densities for d — » oo at T = 0.5 
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FIG. 5. Time series of the ensemble average < i;^ > as a function of time t for different d 
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FIG. 6. Relation between average velocity squared in the frame moving with the current, < > — < Vx >^ , and parametric 
temperature T for the infinite dimensional model. Equipartitioning of energy would imply < > — < Vx 2T. The 

numerical uncertainty of each point is less the size of the symbols 



16 




17 



0.40 \ 



0.36 



0.32 



0.1 



0.5 



0.9 



FIG. 9. Conductivity a{e) as it varies with field strength e for T = 0.5. The curve consists of 90 data points, the numerical 
uncertainty of each point is less the size of the symbols 
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FIG. 10. (a),(b) Poincare section of {/3,sm^) defined in Fig.|l] at the moment of the collision for field strength e — 1. In (a) 
our model with the baker map has been used, in (b) the baker has been replaced by a random number generator, (c) Poincare 
section of l3 at the moment of the collision for varying field strength e.(a).(b),(c) with T — 0.5. 
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